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conformal. We measure the /3-function and the coupling constant as a function of the energy scale. 
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I. INTRODUCTION 



In preparation for the phenomenology at LHC, there 
has recently been increased interest in field theories which 
appear (quasi) conformal when probed towards infrared, 
i.e. over large distance scales. There are basically two 
novel applications which utilize such theories, namely 
walking technicolor [l[ and unparticles 0, . For model 
building along these directions it is very desirable to sin- 
gle out candidate theories which allow for the existence of 
a nontrivial infrared fixed point and then use these as a 
basis for more detailed analyses. To determine the loca- 
tion of the conformal window as a function of the number 
of colors, flavors and fcrmion representations is a non- 
perturbative problem. Traditional way to obtain semi- 
quantitative estimates is the so-called rainbow approxi- 
mation to the Schwinger-Dyson equation for the nonper- 
turbative fermion propagator 0, [1, @ . More recent devel- 
opments include analytic /^-function ansatze and 
analyses on the role of topological excitations ■ Phe- 
nomenological constraints on walking technicolor require 
that conformal behavior is obtained with only modest 
amount of new matter fields, and utilizing higher rep- 
resentation fermions, several plausible candidates have 
been proposed in [ll|, [T^]. Although subject to uncer- 
tainties, these methods are invaluable in providing hints 
towards theories which are of phenomcnological interest 
and motivate their analyses using lattice methods. 

Guided by the phase diagrams of [ll|, initial lattice 
studies of two and three color theories with two fermion 
flavours in the two-index symmetric representation have 
recently appeared [H, [ll, [H, [H, [11. For related 
studies with fundamental fermions in SU(3) gauge theory, 
see [11 in mill [H. 



In this work we study the evolution of the coupling 
constant of the candidate model with "minimal" field 
content, consisting of the standard SU(2) gauge field and 
two massless Dirac fermions in the adjoint representa- 



tion. The euclidean Lagrange density is 
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generators of the adjoint representation. To obtain the 
evolution of the couplig we use the lattice Schrodinger 
functional scheme, which has been very successfully ap- 
|)l^d to e.g. two-flavour QCD by the Alpha collaboration 

Perturbatively, the two-loop /3-function of this theory 
with Nf adjoint fermions is 
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where (3q ^ (22 - 8A^/)/3 and Pi = 4(34 - 32iV/)/3. 
For Nf = 1,2 the coefficient Po is positive while /3i is 
positive for Nf = 1 but negative for Nf = 2. Hence, 
for Nf = 2 there is a zero of the two-loop /3-function 
at gl = -(/3o//3i)16^2 ^ 7 9^ -^ ^ ^ ^2/(4^) ^ q q^ 

corresponding to an infrared fixed point (IRFP). Large 
value of gl casts doubt on the validity of the perturbative 
result and nonperturbative measurement of the running 
coupling on the lattice is required. Indeed, as we report 
below, the lattice simulations reveal large deviations from 
the two-loop perturbative behaviour. 

We note that the /3-function has been calculated up 
to four loops in the Minimal Subtraction (MS) scheme 
[25I I. Beyond two loops the /3-function is scheme depen- 
dent and the result cannot be directly compared with the 
lattice Schrodinger functional scheme used here. 

For the gauge group SU(2) considered here the two- 
index symmetric representation coincides with the ad- 
joint representation. We have recently performed large- 
volume high-statistics lattice simulations of this theory, 
concentra ting on the chiral properties of the mass spec- 
trum [2^, [23]. The main result of this study was that 
at small enough lattice gauge coupling, in contrast to 
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QCD, wc did not observe spontaneous chiral symmetry 
breaking as the quark mass was decreased. Because the 
chiral condensate generates a mass scale, its existence 
would exclude the IRFP. While the non-observation of 
the spontaneous chiral symmetry breaking is very sug- 
gestive, it by no means proves that the IRFP exists at 
non-zero coupling, and one needs to measure the evolu- 
tion of the coupling directly ^. 



II. LATTICE FORMULATION 

The lattice theory is defined with the action 

S = Sq + S'f, (3) 

where Sq is the standard Wilson plaquette action for 
SU(2) gauge fields: 

Sg^(3lJ2 (l - ^Tr[C/,,^C/.+^,,[/]^,^C/t . (4) 

S'f is the Wilson fermion action for two Dirac fermions 
in the adjoint representation of SU(2): 



/=u,d x,y 



(5) 



where 



M^y = S,y-^J2 + Ta')^-.^' + (1 - • (6) 

It differs from the standard Wilson fermion action only by 
the replacement of the fundamental representation link 
matrices Up,{x) with the adjoint representation ones: 



V;:\x) = 2Tr(A'^C/^(x)A''C/t(x)), 



(7) 



where A° = ^"'"i « — li2,3, are the generators of the 
fundamental representation. The lattice action is con- 
ventionally parametrised with two dimensionlcss param- 
eters. 
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where a is the lattice spacing, and and rugfi the bare 
gauge coupling and quark mass. Because the Wilson ac- 
tion explicitly breaks the chiral symmetry of the contin- 
uum quark action, quark masses are additively renor- 
malised and the bare quark mass must be tuned to 
achieve the massless limit. In order to keep the method- 
ology simple, we do not implement 0(a) improvement on 
the lattice. 



In earlier studies [l3|, |26[ it has been found that in this 
lattice theory there exists a critical coupling (3l,c ~ 2 so 
that when /3l < /3l,c the theory shows signs of sponta- 
neous chiral symmetry breaking, but reaching the limit 
niQ is prevented by the appearance of a first order 
transition where mg jumps to negative values. However, 
when (3l > Pl,c the mg — > limit can be reached with- 
out observing the spontaneous chiral symmetry breaking. 
This is the region where it is possible to find conformal 
behaviour in the limit of massless fermions. 

Because the Wilson fermion action breaks the contin- 
uum chiral symmetry, the quark mass is additively renor- 
malised and we determine the quark mass through the 
PCAC relation 

2mQ = [dt{VA{x))/{Vp{x))]^^^^^^^^ (9) 

where Va and Vp are (non-improved) axial and pseu- 
doscalar current correlators [2^: 

VAix) ^a^Yl (d(a:)7075u(x) u(y, 0)75d(z, 0)) (10) 
Vpix) ^a'^Yl (d(a;)75u(x) u(y, 0)75d(2;, 0)) (11) 



Here d and u refer to two quark fiavours, and at time 
slices xo = and xq = L the gauge fields are fixed to 
specific Schrodinger function boundary values described 
in detail below. For each lattice coupling /?l we tune k 
so that TTiQ vanishes. This determines the critical line 
Kc{Pl)', from now on we assume that we always tune to 
the limit mg = 0. 

The coupling is measured using the Schrodinger func- 
tional (SF) method [13, HH, [33, i-e. introducing a 
background field using special boundary conditions and 
studying the response of the system to changes of the 
background field. The scale at which the coupling is mea- 
sured is determined by the finite size of the system. Here 
we consider lattices of volume V = = (Na)^, and fol- 
lowing (33| . the spatial gauge links on the xq = and 
Xq = L boundaries are fixed to constant diagonal SU(2) 
matrices: 

Uf,{xo = 0) = cxp{-ir]a3a/L) (12) 
Uf,{xo^L) = cxp(-i(7r - ?7)cr3a/L), (13) 

where is the third Pauli matrix. The spatial gauge 
field boundary conditions are periodic. During the sim- 
ulation the fermion fields are set to vanish at boundaries 
Xq — and Xq = L, and are "twisted" periodic to spatial 
directions: %jj{x + Lii) — cxp{in/5)il'{x) (3^ . 

At the classical level the gauge field boundary condi- 
tions generate a constant color diagonal chromoclectric 
field, and the derivative of the action with respect to ?/ is 
easily calculable: 



^ Alternatively one can study the properties of the physical degrees 
of freedom in the limit of zero quark mass as outlined in 1281 . 
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mga 


2.05 


0.18625 


-0.00377(29) 


2.2 


0.1805 


0.00016(16) 


2.5 


0.17172 


-0.00079(10) 


3 


0.161636 


0.00213(11) 


3.5 


0.155132 


0.00028(11) 


4.5 


0.14712 


0.00000(05) 


8 


0.136415 


-0.00038(04) 



TABLE I: Values of I3l and the estimates of used. Kc is 
determined on 16* lattices with small statistics runs. The ta- 
ble also shows the small residual value of mga measured from 
the much larger statistics production runs on 16* lattices. 

where, on a lattice of size ~ (iVa)* (ssj . 

fc^-24(L/a)2sin^-^. (15) 

At quantum level the boundary conditions generate a 
background field. Using an effective action F, 

e-r = I I?[C/>>]e-^, (16) 

the coupling constant is defined through 

^ = /^\ = ^ (17) 

We fix rj = 7r/4 after taking the derivative. For our 
action the observable {dSG/di]) is proportional to the 
expectation value of the boundary plaquette and easily 
measurable. The result is the coupling constant g^{L) 
as a function of the physical size of the lattice. Because 
a priori the lattice spacing a is unknown, we denote the 
lattice measurements with g'^{L/a, Pl)- Formal contin- 
uum limit is obtained by keeping g^ constant while taking 
L/a ^ oo. 

III. RESULTS AND ANALYSIS 

As stated in the previous section we will consider lat- 
tice couplings Pl < Pl,c ~ 2 and tune to the uiq = 
limit corresponding to a critical value Kc{(3l)- The val- 
ues of well as the corresponding mq , at each value 
of used in this study are shown in Table |T1 The /Si- 
values which we consider are all larger than the critical 
value I3l,c ~ 2; indeed, already at Pl = 1-9 it was impos- 
sible to reach mg = limit. 

The measurements of g^{L/a, Pl) are done on lattices 
with volumes V/a^ = - 20^, with the results shown 
in Fig. [T] and tabulated in Table IIIII We use a hybrid 
Monte Carlo (HMC) algorithm, with trajectory length 
~ 1; the table also shows the number of HMC trajectories 
for each point. The measurement of the coupling is done 
at the end of each trajectory. The acceptance rate in the 





N 




1/9^ 


2.05 


4 


20000 


0.2504(17) 




6 


51200 


0.2233(16) 




8 


29000 


0.2179(36) 




12 


27000 


0.2187(33) 




16 


30900 


0.2391(53) 




20 


17300 


0.257(22) 


2.2 


4 


16000 


0.3069(17) 




6 


51200 


0.2857(15) 




8 


19900 


0.2719(36) 




12 


55000 


0.2764(30) 




16 


30141 


0.2857(54) 




20 


8000 


0.284(17) 


2.5 


4 


15000 


0.3990(16) 




6 


51200 


0.3815(14) 




8 


21500 


0.3653(32) 




12 


37400 


0.3621(39) 




16 


33000 


0.3618(58) 




20 


4300 


0.369(14) 


3.0 


4 


12000 


0.5427(17) 




6 


51200 


0.5193(14) 




8 


15500 


0.5113(44) 




12 


48500 


0.4958(36) 




16 


19600 


0.4941(74) 




20 


10400 


0.486(23) 


3.5 


4 


15000 


0.6769(16) 




6 


51200 


0.6495(14) 




8 


17800 


0.6367(37) 




12 


57000 


0.6244(34) 




16 


18700 


0.6105(66) 




20 


4590 


0.613(23) 


4.5 


4 


15000 


0.9370(14) 




6 


51200 


0.9067(14) 




8 


14000 


0.8884(41) 




12 


63500 


0.8849(33) 




16 


33741 


0.8671(52) 




20 


8500 


0.878(19) 


8.0 


4 


12000 


1.8221(23) 




6 


51200 


1.7904(13) 




8 


14000 


1.7619(40) 




12 


44100 


1.7492(41) 




16 


12805 


1.7258(90) 




20 


10000 


1.745(20) 



TABLE II: The measured l/g'^{L) and the number of hybrid 
Monte Carlo trajectories used, A'traj., at each lattice size A'' 
and Pl- 



HMC algorithm remains better than 84% for all cases. 
For comparison, we made also small-statistics study of 
the coupling using two flavours of fundamental fermions. 
These results are shown in Fig. [2l 

What can we conclude from these figures? We remind 
that we can directly compare physical lattice sizes only 
at fixed Pl] different Pl correspond to very different 
lattice spacings. In the continuum L is the sole physi- 
cal scale, and g'^{L) must be monotonic (or equivalently, 
i is a function of g^ through dimensional transmuta- 
tion). Thus, the non- monotonic behaviour of data seen 
at Pl = 2.05 and 2.2 is caused by finite lattice spacing 
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FIG. 1: Lattice measurements of g'^{L/a, Pl)- Continuous 
lines show g^{L) integrated from Eq. 
through lattice points at L/a — 12. 
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FIG. 2: Small-statistic measurements of {L/a, (3l) for 2 
flavours of fundamental representation fermions, to be con- 
trasted with the adjoint representation case in Fig. [T] The 
value of g^ at (3l = 2.2, L/a = 16 is 26 ± 7. Continuous lines 
show g^{L) integrated using 2- loop /3-function, constrained to 
go through lattice points at L/a = 16. 



effects, i.e. too small L/a. Hence, we do not include 
lattices with L/a < 8 in subsequent analysis. 

Our most significant result is readily evident from 
Fig. [H g'^{L/a, /3l) is an increasing function of L at 
Pl Z 3, (g^ < 2.2), but decreasing at Pl ^ 2.5 (g^ > 2.2) 
at large volumes. This is precisely the expected be- 
haviour for a theory with an infrared fixed point and 
markedly different from the asymptotically free QCD- 
likc behaviour in Fig. [51 In this case g^{L) is always 
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FIG. 3: Step scaling A{Li, L2, Pl) for different pairs of vol- 
umes, see Eq. p8|l . plotted against 1/ g^ {L\, (3l)- The scaling 
factor is 2, except for the pair with largest volumes. Also 
shown is the 2-loop perturbative result for corresponding ra- 
tios of volumes. 



increasing, and the larger is, the steeper the increase. 
For the adjoint representation the observed fixed point 
2 . . . 3 is significantly below the two-loop result 7.9. 
Further analysis is needed in order to estimate the full 
/3-function. The standard method proceeds through the 
measurement of step scaling functions: using lattices of 
sizes L and sL, where s is a scale factor, we measure 
g'^{L/a, (3l) and g'^{sL/a, (3^), thus obtaining a discrete 
step of the evolution of g^{L). Now we can search for 
P'^ so that g'^{L/a, f3'^) = g'^{sL / a, Pl)- Repeating these 
scaling and matching steps the evolution of g^ over a very 
large range of scales can be covered. The continuum limit 
extrapolation is enabled by redoing the analysis on differ- 
ent original L/a [11]. This method has been successfully 
used, e.g. to determine the running coupling in QCD 

In Fig. [3] we show the step scaling using the quantity 



A(Li,L2,/3l) 



1 



1 



g2(i2,/?L) g^{L,,pLy 



(18) 



We use scaling factor s = L2/L1 ~ 2, except for the 
largest volume L/a = 20, where we do not have L/a= 10 
data avalable and use s = 5/3. We also show the 2-loop 
scaling result computed from Eq. 

If finite lattice spacing effects were absent, all points 
with the same s should fall on a universal curve. This 
is clearly not the case; the small volume data lies signifi- 
cantly below the large volume points. We do not attempt 
to extrapolate the data to infinite volume (continuum 
limit), because we do not trust the smallest volume data 
to be close enough of continuum behaviour. Nevertheless, 
A(Li, L2, Pl) becomes negative at small l/g", indicating 
the existence of the infrared fixed point. 
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The main reason for the difRcuhy of using the step 
scaUng is that the evolution of g^{L) is very slow, much 
slower than with fundamental fermions. The slow evolu- 
tion is easily masked by finite lattice spacing artifacts at 
small i/a, whereas the faster evolution in fundamental 
fermion theory is rather easily observable. Improved ac- 
tions can be expected to help significantly here. The slow 
evolution also increases the number of required recursive 
steps in the step scaling method and hence statistical 
errors. 

Due to the numerical problems with the step scaling 
function mentioned above, we will therefore also con- 
sider the following alternative method: we use the mea- 
sured values of g^{L/a, Pl) on volumes 12'*-20^, and fit 
a (3-function ansatz to the data. The ansatz includes 
the small-f/^ perturbative part and it should have only a 
few tunable parameters. The aim is to find an envelope 
of continuum /3-functions which arc consistent with the 
measurements. The underlying assumption is that the 
lattice artifacts are already sufficiently small at i/a = 12. 
This should be checked with new simulations with im- 
proved actions in the future. 

The simplest ansatz is the perturbatively motivated 
one P{g) = —bog^ — big^ — 629 ''• Here bo and 61 are fixed 
by the perturbative /3-function Eq. ([2]), and 62 is a fit 
parameter. This function is actually consistent with the 
data, but it is very likely too constrained giving unreal- 
istically small error range. In order to have independent 
parameters for the value of the fixed point coupling and 
the slope of the /3-function at the fixed point, we are led 
to the following ansatz: 



0.08 



0.04 



-0.04 




FIG. 4: The /3-function obtained from Eq. HHJ, with = 2.2. 
The shaded area shows the estimated error range of the fit. 
Shown are also universal perturbative one- and two-loop /3- 
functions, together with the three- and four-loop results in the 
MS-scheme Q^. Because of the different scheme, these are 
not directly comparable with the lattice (SF-scheme) results. 



(.9*>'5) ~ (2.0,6) to (3.2,15) in a narrow region. Thus, 
the parameters are strongly correlated. The slope of the 
/3-function at the fixed point is a universal quantity, but 
it is only weakly constrained: 



Pig) = -bog' - big"- + {bogt' + b,gt-')g 



(19) 



Again &o a-nd bi are fixed by the perturbation theory and 
g^ and 5 are fit parameters which determine the fixed 
point and the slope of the /3-function around the fixed 
point. The ansatz is naturally not of asymptotic form 
at small 5^, but this docs not matter for our purposes: 
the g^-tcrm turns out to be negligible at small g^, and 
the purpose of the ansatz is to approximate the true (3- 
function in the vicinity of the fixed point. In practice 
equivalent results are obtained using perturbatively mo- 
tivated ansatz ...625^ + &3ff^; however, Eq. offers 
more freedom for the slope at g'^ = g^. 

We perform the fit as follows: integrating (3{g) = 
—Ldg{L)/dL and constraining g^{L) to the lattice mea- 
surements at L/a = 12 for each Pl gives us the contin- 
uous lines shown in Fig. [if. The integrated g'^{L/a, Pl) 
are compared with the lattice measurements at L/a ~ 16 
and 20, giving us a criterion for the fit. Good fits 
(x^/d.o.f~ 1) are obtained when parameters vary from 



^ Strictly speaking the fit is not optimal due to the constraint at 
L/a = 12. Relaxing this constraint gives better x^t cost 
of a more complicated fit procedure. However, because L/a = 12 
has smaller errors than larger volumes the effect is small. 



' dg 



0.12... 0.42. 



(20) 



=9i 



We take 5^ = 2.2, 5 = 7 to be our benchmark value, 
which is motivated by the fact that S ~ 7 is the or- 
der of the NNLO perturbative /3-function. These lines 
are shown in Fig. [T] However, we by no means imply 
that the fitted value gives the perturbative result; it only 
quantifies all unknown contributions. If we fix (5 = 7, 
can vary between 2.05 and 2.4. 

The result of the fitting procedure is shown in Fig. 21 
together with the error band formed by the variation of 
the parameters in the region where acceptable fit is ob- 
tained. While the error band is wide, the qualitative 
features agree with the corresponding step scaling func- 
tions calculated from the measurements; both methods 
indicate positive f3{g) at large g^. We also note that the 
two-loop perturbative /3-function alone fits the measure- 
ments weU at /3l > 3 {g^ < 2). 

Fig. H] al so shows the three- and four-loop MS-scheme 
results [2g. Because of the different scheme, these are 
not directly comparable with the SF-scheme results ob- 
tained in our simulations — for example, the value of g^ is 
scheme dependent. However, it is interesting that also in 
the MS-scheme the fixed point is at substantially smaller 
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FIG. 5: g^ifj. = 1/L) determined from Eq. (HSJ, together 
with perturbative g^. Scale Ao is determined so that ^ — Ao 
is at the IR Landau pole of the asymptotically free one-loop 
coupling and at the UV Landau pole of the measured non- 
asymptotically free branch {g^ > g^). For illustration, the 
L/a = 16 lattice points have been placed on the benchmark 
curve at the corresponding gr'^-values, indicating the relative 
scale hierarchy between simulations. 



than indicated by the universal two-loop calculation.'^ 
Finally, the fitted /3-function can be integrated to ob- 
tain the coupling g^. This is shown in Fig. [5] as a func- 
tion of /i = 1/L. The result is split into two branches: 
the asymptotically free one with g^ < gl, and the not 
asymptotically free with g^ > g^. These branches are 
disconnected: when the /3-function is integrated to ob- 
tain g^ (fi) , the integration constants are independent on 
both sides of the fixed point. Thus, the scaling between 
the /i- values between the branches is arbitrary in Fig. [5] 
Finally, we note that the two-loop perturbative /?- 
function alone gives us an acceptable fit if we ignore the 
data at (3l < 2.5 {g^ > 2.7). Indeed, in that case it is 
possible to fit a /3-function which does not feature an IR 
fixed point. Thus, the evidence for the fixed point comes 
from the data at (3^ = 2.05 and 2.2. These are at rel- 
atively strong bare lattice coupling, and one can expect 
large finite a effects. This can be clearly seen at small 
L/a in Fig.[TJ Nevertheless, as L/a ^ oo at fixed Pl, the 
scale hierarchy between the physically interesting scale L 
and the lattice spacing a increases, and the evolution of 



the coupling constant approaches the continuum one. 



IV. CONCLUSIONS AND OUTLOOK 

Let us now discuss the interpretation of these results. 
We note that the observed range for the critical coupling, 
gl ~ 2.0 - 3.2 is significantly below the perturbative es- 
timate ~ 8 from Eq. ([2]). Similar behaviour has been 
observed in SU(3) gauge theory with 2-index symmetric 
fermions [T^ . The obvious question is why the perturba- 
tive two-loop /3-function fails although the nonperturba- 
tive analysis implies that a, = g^/47r is small. Since the 
perturbative determination of the fixed point is based on 
the competing effects from one- and two-loop orders, it 
is reasonable to expect that higher loop orders can con- 
tribute with similar magnitude and one needs nonpcr- 
turbative analysis. This is seen in the MS-scheme results 
in Fig. m Beyond two loops the value of g^ is scheme 
dependent; our result is in SF scheme. Comparison of 
the nonperturbative result and the perturbative curves 
in Fig.[3] demonstrates that even if the actual value of 
the coupling at the fixed point is small, the nonpertur- 
bative effects can be large. 

When we combine the results of the analysis here with 
the spectrum measurements of ref. [2^ , a consistent pic- 
ture emerges. In the whole parameter range where sim- 
ulations with massless fermions are possible the infrared 
behaviour is controlled by the conformally invariant fixed 
point and there is no chiral symmetry breaking. If 
g^ < g^, the theory is asymptotically free, otherwise not. 
The conformal invariance can be explicitly broken by e.g. 
adding a small fermion mass, which can lead to a "walk- 
ing" type evolution. 

As discussed above, for the lattice action used here 
it has been observed that it is not possible to take the 
limit niQ when Pl ^ 2. Thus, the Pl = 2.05 - 
points in Fig. [5] are close to the smallest Pl where the 
conformal behaviour can be studied. This limitation is 
a non-universal lattice artifact, and may differ for some 
other fermion discretisation. This motivates a study with 
fully 0{a) improved lattice formalism. 
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